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1. INTRODUCTION 

To manufacture various devices of solid-state electronics, heterostructures of different 
configurations are frequently used. For their growth different methods could be used: gas-phase and liquid- 
phase epitaxy, sputtering of materials in magnetrons, and molecular beam epitaxy. The manufacturing and 
use of heterostructures in different devices were described in a large number of experimental works [1]-[11]. 
At the same time, fewer works are describing the prediction of epitaxy processes [12], [13]. In this paper, we 
consider a vertical reactor for gas phase epitaxy as presented in Figure 1. The reactor consists of an external 
casing, a substrate holder with a substrate, and a spiral around the casing in the area of the substrate to 
generate induction heating to activate chemical reactions of decay of reagents and to the growth of the 
epitaxial layer. 

A gaseous mixture of reagents together with gas-carrier inputs into the inlet of the reaction chamber. 
At the first stage of the growth of the heterostructure, a buffer layer was grown on a substrate. Next the 
obtained structure was annealed. After that an epitaxial layer was grown on the buffer layer. The main aim of 
this paper was analyzed of changing of the properties of the final heterostructure under influence of the 
considered annealing after grown of the buffer layer. Additionally, we introduce an analytical approach for 
analysis of mass and heat transfer in multilayer structures with account their nonlinearity, changes of 
parameters of processes in space and time, mismatch-induced stresses. 
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Keeper of substrate 


z=-L ; z=0 z-L 
External casing 


Figure 1. A reactor for the gas phase epitaxy in neighborhoods of the reaction zone 


2. METHOD OF SOLUTION 

The aim of the present paper will be solved by analyzing the Spatio-temporal distribution of the 
distribution of temperature and the concentration of the deposited material. The required temperature 
distribution was calculated by solving the following boundary value problem [14]. 


caen. = p(r,d, z t) + div(A: grad[T(r,$, z €)]-[(, p, z, €) — 0(r, p, z, t|  e(D)- 
T(r, ,z,t)-C(r, ¢,z, t)} (1) 


Where the vector ? is the speed of flow of the considered mixture of gases; parameter c is the capacity of 
heat; function T (7,@,z,2) is the distribution of temperature in space and time; function p (r,@,z,t) describes the 
power density, which stands out in the considered system substrate - keeper; r, 9 and z are the cylindrical 
coordinates; f is the current time; function C (r,@,z,f) describes the distribution of concentration of mixture of 
gases in space and time; parameter / describes the conductivity of heat. The value of heat conductivity could 
be determined by the relation: A = vlc,p/3, where ? is the modulus of the mean squared speed of the gas 
molecules, which equal to ? = ./2kT/m, | is the average free path of gas molecules between collisions, c, is 
the heat capacity at constant volume, p is the density of the gas. 

To solve this boundary problem, one shall take into account the moving of a mixture of gases and 
the concentration of this mixture. We calculate the required values by solving the equation of Navier-Stokes 
and the second Fourier law. We also assume that the radius of the keeper of substrate R is essentially larger 
than the thickness of diffusion and near-boundary layers. We also assume that stream of gas is laminar. In 
this situation, the appropriate equations could be written as (2) and (3). 


T GP - -v (2) + vá, (2) 
200420 — div{D - grad{C(r, $,2,t)] - [BC $,2,t) - Ë]: Cr. $,2,0} © 


Where D is the diffusion coefficient of a mixture of gases (gases-reagents and gas-carrier); P is the pressure; 
vis the kinematic viscosity. Let us consider the regime of the limiting flow, when all forthcoming to the disk 
molecules of depositing material are deposing on the considered substrate, flow is homogenous and one- 
dimension. In this case, initial and boundary conditions could be written as (4). 


C (r.9,-L.t) = Co, C (r.9,0,t) = 0, C (7,0,z,2) = C (7,27,z,0), C (7,9,2,0) = Cod (z+), 


8C(r, S z,t) 8C(r,$,zt) 8C(r,$ zt) 
C 0, „Z,Í * o, = 0, = , T r, X40 = T,, 
(0.9.2.0) oe lees 0$ eco M luo (r,9.2,0) 
aT (7, 9,zZ,t) 4 OT (r,,Z,t) OT (r,9,Z,t) 
Se = oT (R,Q,z,t), ——— EET EE , 
or T=R ( $ ) ag $-0 ag $-2n 
OT(r,$,z,t) dvr (r,S;z,t) 


= = 4 = = 
A 2. bust oT*(r, $, —L,t), T (r.O,zt) = T (r,275z,0), 2 <= 


eee] = mei) etie] E mue] 
ag g=0 09 ag g=0 09 


(4) 
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Portz] o, Q9TO9z0| —— 64 (r, i z, t), T (0,95, + ©, v, (1,0,-L,f) = 0, 
or T=R az z-L 


vr (7,9,0,t) = 0, v (r, d, L,t) = 0, v(r,0,z,t) = vi(r,275z,D, vr (0,9 z,t) £ oo, vo (7,9,0,1) = er, 
Vo (1,9,-L,t) = 0, vo (r,Q,L,t) = 0, vol7,0,2,0 = ve(r,275z,t), V0, qz, t) # oo, v(r.Q,-L,t) = Vo, 
v,(r, 9, 0,t) = v,, v(r.g,L,t) = 0, v(r,0,z,t) = vr,27,z,t), v{0,9,z,t) # ©, v(r,q,2,0) = 0, 
ver, Q,z,0) = 0, v-(7,@,-L,0) = Vo. 


Here parameter c is equal to o =5.67.10° W.m?.K“, parameter T, describes the room temperature, 
parameter @ describes the frequency of rotation of the substrate. Equations for components of velocity of 
flow with account cylindrical system of coordinate could be written as (5a)-(5c). 


Ov, Ov. Vg Ov Ov, 0? v, 0?v. 0?v 0?v. ð (P 

r= y, r_e $ y, z ( Tr r T z) - ( ) (5a) 
ôt or r 0$ Oz ðr? ðrðz = az? oraz or Né 
I ys Ole T NM: 18yr | 2 0jvg — 1 dy, | Ü?vg| i18 G) (5b) 
a Tər r 0$ Z az rürüo r2 ð$?  r?80Q0z Oz? r ab NE 
Ov, Ov, Vg Ov dv 0? v, 0?v. 1 0?v 0 (P 
BE ERE L AME d CIS g, ( NN; NES 2) 2 (2) (5c) 
ôt or r 0$ Oz dz? or? r? 0$? Oz Né 


We determine the solution of this system of equations by using the method of averaging function 
corrections [15]-[20]. Framework this approach to determine the first-order approximation of components of 
the speed of flow of a mixture of gases we replace the required functions on their average values v,>a1,, 
Vo Qip, V 05; in right sides of equations of system (5). After the replacement, we obtain equations for the 
first-order approximations of the components. 


= a(R) a= Oe = OD) n 


Solutions to the above equations could be written as (7). 


ð rtP 
Lb Fat, Vig = dirt v=- hdr (7) 


The second-order approximations of components of the speed of flow could be obtained by 
replacement of the required functions on the following sums v,>a@1,, Vg> 019, V Qız. Approximations for 
the components could be written as (8a)-(8c). 


avor _ RM avr var oe) a aye rr (azp *v19) vir 
at dp Tae oz? | Droz ar VE (a, + vir) r r ag 
Ov. 
“(Ges + viz) z 2 (8a) 
0v2$ 18?v,. , 2 0?v,g — 1 O0?v,, , Vig 10 (P a 
= T = — (a5, tv Vir) 
ôt rürüó  r? Ag? r? 0$0z 0z? rôp 
(a2g+19) vig ut 
CÓ T ME — (a, + Viz) (8b) 
Ova, _ Re z1 v) ð 2 (2) - (an +v ) ee (a2ptvig) viz 
at az? ðr? — r2 ag? az V 25 ACTAE T m 
Ov. 
—(@2z + viz) S (8c) 
Integration of the above equations leads to the following result. 
- a? vir 3vir vir A) ð ( tP ) t ðvir 
=v, ( ar2 t oróz az? T ürüz dr ar Jo H dr Jo (a2, + Var) ar dt 
t (@2gtV19) 0vir t Ovir 
zi ET oe IT Sq ae + Viz) 57 dt. (8d) 
24e 255a, ire Pep P Gs) 
T ae M 0$? r? 0$0z Oz? rôp 
t (a +v jm 
T (@zr sk vr)? at dt -Í eet dt -f (a2; T viz) t Se dt, (8e) 
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a? viro 0?viz 1 0?v4z ð tP t ðviz 
= -vf e a ta FE Jar - Z (f zdr) — fa oe + VAS) 2s diet 


t + üviz t avaz 
piema a, + vio) 55 —t dr- Í (@2, + viz) x dt (8f) 


We determine average values o»; Q29, 00; by (9). 


dS fo So So” faa, — vidzdodrat, 


TOR?L 
1 O çR 2 L 
tog = mh hr h SiOz- vis)dzdodrat, 


fo So K E Oz- Miz)dzdgpdrdt, (9) 


a — 
2z sum 


Where © is the continuance of moving of the mixture of gases through the reactor? Substitution of the first 
and second-order approximations of the required components of speed into the relation (9) give us the 
possibility to obtain a system of equations to determine required average values. 


A105, t B4026 t C40; = Di 
4205, + B505$ + C505; = D; (10) 
A305, + B3026 + C3057 = D3 


Where A4, = 1+ fj(6  O fo rf, fi," dzdódrdt, B, = fj (@ -O fo fo" f, er dz dódrdt, 


= = 2p2 _ 2n Om, vir 3vir , O0?vaz 
éco usu v pr E( r p ar Paru ded dr — 


[2] R 2n cL O FR r2n pL 
T 92 p22 OV. QU, 
—t)dt — g? RV -| (0-0]r Us gc dnnt — ne ag dzdd¢dr x 
0 


x (0 — dt, A; = f°(@— 2) frr pen i ? d; dá drdt, B,-14f? p jy JE 9n d zdġdr(0 — 


2 10?y, , 2 0?v 10?v,, , 0?v 
-t dt, D, = v ff ape ere p mm +2 RE urat =) Dnm 
-fe-O0ferfr [vie m dzdodrdt -20?R?V? — fj (0 -O fy fo" feu vs mar dad drdt, 
9 a 


As = f (0 — t) p E. o. 2 dzdpdrdt, B, = f (0 -DK rf, f. Pou panier C;=1+ 


L 0o 
+= 6?R?Vp, D, = v Jf (8 — of nde (Er + Pu a i a) dzdpdrdt — ff (© — t) x 
Mi uf [os 77: dzd$drdt — f; (0-0f rf vs 7 dzdódrdt — 10? RN. 


The solution of the above system of equations could be deemed by standard approaches [21] and 
could be written as (11). 


Qr = A,/A, Q6 = A4/A, Q2; = 4z/4, (11) 


Where A = A, (B2C3 — B3C2) — By(AzC3 — A3C2) + C,(A2B3 — A3B2), A, = Dy (B2C3 — B3C2) — 
—B, (D2C3 — D3C2) + C, (D;Bs — D3Bz), Ag = Dı (B2C3 — B3C2) — Bı (D2C3 — D4C2) + Cy X 
X (D2B3 — D3Bz), Az = A1(B5Ds — B3D2) — Bı (A2D3 — A3D2) + Di (A; Bs — A3B2). 

In this section, we obtained components of velocity of the stream of a mixture of materials in the gas 
phase, which are used for the growth of heterostructure, and gas-carrier in the second-order approximation 
framework method of averaging function corrections. Usually, the second-order approximation is enough 
good approximation to make a qualitative analysis of obtained solution and to obtain some quantitative 
results. Let us rewrite (1) and (3) by using a cylindrical system of coordinates. 


regan ES [Pree 4 E D + regan pui 2 (c(r, Pzt) 
T(r, ¢,z,t)-[v,(r, 6, z,t) — č, (r, ġ,z,t)] j- £35 [vor $.z, PLAN $,2z,t)|-C(r,¢,z,t)- 
T(r, ¢,z,t)}—c- zv "«(rnd;zu-vmdq0üzt]-C(rózt-TGmróz0l)-p(rézt (12) 
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OC(7,¢,z,t) = 18 [ro akero) + 


E 12 ds “|p acea] + Z [p eee} - 2 ? (rC(r, $, z, t) - 
[v,(r, b, z, t) - v, (r, d. z t)]} — - rcr, $,zt) [vo (r,d,z, t) — (r, ġ,z, t)]} — 
2- (Cr, d, z, t) a [v,(r, $, Z, t) ~ D, (T, $, Z, t)]}. (13) 


ror 


We calculate the distribution of temperature in space and time and the same distribution of 
concentration of the mixture of gases we used the method of an average of function corrections. To determine 
the first-order approximations of the required functions we replace them on their not yet known average 
values ar and aic on the right sides of the above equations. Further, we used a recently consider algorithm 
to obtain the first-order approximations of the temperature and concentration of the gas mixture. 


TG, $, z, t)= = T, iE [MA dt — a, rec IN E e 


E tovg up (22) t 0[vz(r,G. zt) -vz(r.duz,c)] 
1T oh $ 5 $ dT — G4704c tf dr (14) 
exc (t8trlvr(rdizm) vy .b.z,7)]} arc pt êlor pz) -ier pz] 
r ane A a eae is 
The above not yet known average values could be calculated by using the standard relations. 
R 2m pL 
= er Sale Ío r f, Lm T,(r, $, z,t)dzdddrdt, 
L 
= — N i f ie Ci(r, b, z, tT)dzdgdrat. (16) 


Substitution of the first-order approximations of temperature and concentration of gas mixture into relations 
(16) gives us the following results [20]. 


aic = Co/L- {1+—— J (0 — t) f?" f^ [v.(R, dz, €) - 9, (8, dz, t)]dzdpat 4-53], 
ar =|T, +f O -DN p JO" JE Bez ay d$drat] (1 - ——— fete - 0 x 
x JÈ" ^ In (dz. — 8, (Re daz, Dldzdpdt — [e JË f?" flo. (r 4.2.2) - 

2, (r, d, z, t)]dzdġdr x (0 — t)dt — + vo} = fio - 0f flv B20) - 
b,(R, p, zt) |dzdpdt + 1+ GV,/RLY-?) 


The second-order approximations of temperature and concentration of the mixture of gases 
we determine the framework for the method of averaging of function corrections [15]-[20], i.e. by 
replacement of the required functions on the right sides of (12) and (13) on the following sums 7 oor +T), 


C— ac +C\. In this case, the second-order approximations of the above-required functions could be written 
as (17a) and (17b). 


c: To (ds z, t) =a fi” HOWE egi A us need gr seu Hune 

=g IKI v, (r, p, z, T) — v.(r,d, z1)]: [aoc + (r, ġ,z,T)] : [aor + T4 , z,t)]}dt 

z 5 fo {lec t GG 4,22): [vs(r 4,2, — Bg (7, p, zT] - [eor + T(r, di z, ])dt 

-c-2 flog, dz, — 9,0, 2,0] [aac + C(t do] - lazr +T, Gs dz, ]de 

+f pr, d. z, dt T, (17a) 


Co(r, dz, t) = Lr rp 2500929 dr TE als DBE a tt+of , D O20 ar 

12 fr flee +0820 fl - n hou] 

om Lh læzc + Clr, p, z,2)] : [vo (m p, 2,7) — v3 7, z, )]dv + Cod(z + L) 

ES IN [a2¢ + C, $.z, 1)] i [v Gr, $, Z, T) p v,G, $, Z, Tt)]dī. (17b) 
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Averages values of the second-order approximations of temperature and concentration of mixture 
oor and aoc have been calculated by using the following standard relations. 


2n 
Gu = —zh h rh [T TOdzdedrat, 


0 cR 2 L 
ax = —z h fork, SC — COdzdédrdt. (18) 


Substitution of the first- and the second-order approximations of temperature and concentration of mixture 
into relations (18) gives us the possibility to obtain equations to determine required average values. 


damian k — D" TE 6,2, Ddrdodt — [j(6 — 
t) fi A, TQ, bz, t)dzdġdt x— LL tL (0 = t) f; f, T. (0, 6,2, t)dzdódt — 

p (9 —t) ZÉ ATR, Z, a x [azc + 6,05, 6,7, t)] — es 70a 0v (R, , 2, t) — 

(R, $, Z, t)]dzd$dt —— — — xf i rf ,Ui(r oz t)læz + Clr, o, z t)] - arac}: 
[v-(r, b,z,t) — 9, (r, p, z, t)]dzdġdr x (0 — t)dt — Eh (0 — of F^ [TA (r, by L, t) (ac + 
Co) — 4704 :]dódrdt = xf (0 -t) jt [v,(R, $, z, t) — V(R, Q, z, t)]laac + 

C,(R, d, z, t)|dzdpdt + 1-—_ x fr (-o0f rf Er, gz t)’, pz t): 
[aoc + C4 (r, ġ,z, t)ļdzdġdrdt + (az¢ + Co) x +20V,/LY +, 


2 OC, (7,9,z,7) 8C4(r, zx) e 
tre = urhe (8 — OI Sy” DE a e alpara foco - 


t) IW „{rlazc — aic + C (0,2, 0)] : [o Q5 6, z, 0) — 9 (R. d z,t)|}dzdodt x —— 
j- (0 — t) f; A "(ce — ic + Cy) dzdgdrat. 


— TE 
T 
After growing the buffer layer, we consider annealing the resulting two-layer structure. During the 
annealing, one can find diffusion mixing of the heterostructure layers. Thermal diffusion in this case will be 
one-dimensional and perpendicular to the interface between the layers. But due to the mismatch of the lattice 


constants of these layers, one can find mismatch-induced stress. In this situation to describe the mixing of 
layers, we use the second Fick's law in the following form [22], [23]. 


ap(xyzt) | 0 T dp(xyzt) Lz 
S = 2 |p ell +2— 2 |?s x Vsulx,y,z, t) J, pGo y, W, eaw] + 


ð 
07 ES vsu yz, o "S taw] ae 


With initial and boundary conditions. 


Op(x,y,z,t) =0 Op(x,y,z,t) ET ee zu Op(x,y,z,t) =0 
ax x=0 f ax x=Ly > i ay x=Ly , 
Op(x,y,zZ,t) Op(x,y,zZ,t) 
res = 0, ——— = 0, p (wy,2,0)=fp (x,y,z) 
Oz z=0 Oz 


x=Lz 


In the above relations, the following denotations were introduced: Q is the atomic volume of the 
MAE : i i ‘ . hz : 
dopant; V; is the symbol of surficial gradient; o (x,y,z,T) is the diffusant concentration; if; p(x, y,z,t)dz is 


the surficial concentration of the considered diffusant on the interface between layers of heterostructure; 4 
(x,y,z,f) is the chemical potential due to the presence of mismatch-induced stress; D and Ds are the 
coefficients of volumetric and surficial diffusions. Values of dopant diffusions coefficients depend on 
properties of materials of the heterostructure, speed of heating and cooling of materials during annealing, and 
Spatio-temporal distribution of concentration of dopant. Dependences of dopant diffusions coefficients on 
parameters could be approximated by the following relations [23]. 


YGoyz,t) YGoyzt) 
De = Di (x,y,z, T) [1 + § PO], Ds = Du Gv D [1 + 6 SPO (20) 
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Here Di (x,y,z,T) and Dzs (x,y,z,T) are the spatial (due to accounting for all layers of heterostructure) 
and temperature (due to Arrhenius law) dependences of dopant diffusion coefficients; T is the temperature of 
annealing; P (x,y,z, T) is the limit of solubility of dopant; parameter y depends on properties of materials and 
could be an integer in the following interval y e[1], [4], [23]. The concentration dependence of diffusion 
coefficients was described in detail in [23]. Chemical potential 4 in (19) could be determined by the 
following relation [24]. 


4 -E(z)Qoi; uiGxy zt) euo y cst) V2, (21) 


* E ð x s 
Where F(z) is the Young modulus, oj is the stress tensor; uj; = =3(# + a) is the deformation tensor; ui, uj 
J Li 
are the components ux(x,y,z,f), Uy(x,y,z,t) and u-(x,y,z,f) of the displacement vector U(x, y, z, t); xi, x; are the 


coordinate x, y, z. The (21) could be transformed into the following form. 


* 2|89ujGyzt) , 9ujGoy 0 | [1|89ujGyt) , 9wjGoyz0| — 
i, 0 y; z t) D E(z) 2 | ax; + Ox; 2 Ox; i Ox; o9; 


o(2)óij [QukGoy,z.t) 
d [RZAD - 3e] - KEBO Œy, z t) - To] dij}. 


. [|8ujGoyzt) , 9ujGoyz.t) 1 |9uiGoyz,t) du; Gy, z,t) 
ux, pt z t) F | Ox; + OX; 2 + 


3eo] - K@)B@MTy, z,t) - Tolôy} 7 EO). 


o (2)óij [2 (x,y,Z,t) m 
1-2o(z) 


Ox; Ox; Oxk 


Where ois Poisson coefficient; & = (as-agr)/agr is the mismatch parameter; as, az; are lattice distances of the 
substrate and the epitaxial layer; K is the modulus of uniform compression; fj is the coefficient of thermal 
expansion; 7. is the equilibrium temperature, which coincides (for our case) with room temperature. 
Components of displacement vector could be obtained by solution of the following [24]. 


a? Ux(%,Y,Z,t) _ oxx(%,y,Zz,t) , 00xyGGY,Zt) , Aoxz(x,y,z,t) 

pn—u 0t? = ax F ay T az 
a? uy(GXyzt)  O0eyxGuyzt) | Adyy(%yz,t) , 00yzGoy,z.t) 

p(z) at2 ~ ax T ay + az 
07 uz(x,y,z,t) _ 00zx(x,y,Z,t) dozy(x,y,Z,t) 00zz(x,y,z,t) 

PO Ga at2 ax T ay E az 

Where 
_ EZ) duj(x,y,z,t) , Ouj(x,y,z,t) " óij Quy CGy,z,t) 3 Quy (x,y,z,t) 
oy = es [uen 4 URED _ Su Penne + cra, x MEP BOOKET, y, z t) — Fy] 


p (2) is the density of materials of the heterostructure, ó; Is the Kronecker symbol. With account the relation 
for oj; last system of the equation could be written as (22). 


2 9? 
pz uy Ux(x,y,zZ,t) se )4 5E(z) }2 Ux (x,y,Z,t) 4 (kc || E(2 } x Uy (x, y,z,t) 4 


0t? 6[1+0(z)] 0x? 3[1+90(z)] axdy 
E(z) [^ uyCGoy,zt) O0?uz (x,y,z,t) | E(z) l 07uz(x,y,z,t) x. T y,zt) 
iesu. y T az KO) rrt ae OR ee 
a? E yzt) E(z) [Len 07 uy (x, y,z,t) - ?re yzt) 
p ) ^ 2[1+0(z)] Ox? T axdy ^ K(z)B (z) + 
à { m jen + duz(x,y,z,t) } 4 07 uy(x,y,2,t) 
z (2[1+0(z)] Oz 0y? 
5E(z) } ( || E(2 ) 07 uy(x,y,z,t) 0?uyGoy,z,t) 
% A + K(z) + K(z) 6[1+0(z)] Oydz T K(z) Oxdy 
( uen = E(z) [eee 07uz(x,y,z,t) 07 uy (x,y,z,t) 0?uy(x,y,z,t) 
P at2 ^ 2[1+0(z)] ax2 8y? axdz ayaz 
Oux(x,y,z,t) , 0uyGoy,zt) DU YZ, 2 ++ = { E(z) | Ouz(x,y,z,t) e Oux(x,y,z,t) = 
Oz KG ) | Ox F ðy T 6 0z \1+0(z) Oz Ox 
Ouy(Xy,zt) duz(x,y,z,t) arya 
omm Mazan) kopa reae (22) 


Conditions for the system of (8) could be written in the form. 
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0u(0,y,z,t) Ou(Lx,y,z,t) Ou(x,0,z,t) Ou(x,Lyz.t) 
= 0; = 0; = 0; = 0; 
Ox Ox oy oy 
du(x,y,0,t) du(x,y,Lz,t) 22 END E 
a; Oz m 9 UY, 2,0) = Uo; U(X, y,z, oo) = to 


The distribution of the concentration of the diffusing substance has been calculated by using the 
method of averaging functional corrections. The first-order approximation could be determined by the 
following relation. 


a rt Esa! 
pi(%, y, z t) = Mc A= ho Ds, Go y, z, T)Vsu (x,y,z, Déli + a dt + 


ac — in Ds, Go y, z, T Vsus Qe, y, z, v) É zhi + | dt + f, (x, y, Z) (19a) 


PY (x,y,z,T) 


The average value of the considered approximation of the considered function could be calculated by using 
the standard relation [15]. 


1p = —— fr f? es je p1(x,y, z, t)dzdydxdt (23) 


zh 


Substitution of ielapon (19a) into relation (23) allows obtaining the desired average values in the 
f zi Lr ty iA z foGoy,z)dzdydx. Next, we calculate the second-order 


following form: « 
g 1p — erat 


approximation of the considered concentration of the diffusant by using the standard iterative procedure of 
the method of averaging functional corrections [16]. The required approximation was calculated by the 
following relation. 


p2(x,y,z,t) = Z f Dix, yz r)ean f + NE C D + fo (x,y,z) + 


ag 75 V.u(x, y, z, v) f, [azp + p1(x,y, W, c) ]dW dt + a Slice 7S Vou(x, y, z, v) f," [aa + 
kT kT 
P» a y,W,t)|adW dr. (19b) 


The average value of the second-order approximation of the required concentration o», is determined using 
the standard relation [15]. 


azp = —— ff fe Klp, y,z t) — p (x,y,z, t)ldzdydxdt (24) 


Ta 

Substitution of relations (19a) and (19b) into relation (24) gives a possibility to obtain a relation for 
the required average value: 055-0. Next, let us solve equations of system (22), i.e., to obtain components of 
the displacement vector. Equations for the first-order approximations of the considered components after 
appropriate substitution into the (22) take the form. 


p(z y Urx(%yZ,t) _ -K(z Bz ) ened) 


0t? 
p(z) a? MY t) —K(z)B (2) Te t) (22a) 
z 
p(z) 2-9:6»29 — _K (2) B(z 20) 


0t? 


Integration of the left and right sides of the (1b), (3b), and (5b) on time gives us the possibility to obtain 
relations for the above approximation in the final form. 


(z) 8 0 (z) 8 9 
wy 05,2, t) = uos + KG) ESS fo [o Ty, z, ded — KZ) EOS fe f» TG y, z, drdà, 


B(z 8 rt? B(z) 8 po pó 
u4y x, y, z, t) = uoy + K(z) cords J, TG, y, z, 1)dtdd — K(z) Ferre JS, Tay, 2,1) dtd, 
BG 0 BG) a 


w(x, y, Zt) = uoz + K(2) 5 2 f f? TEx, y, z, c)drd8 — K(z) P9 f? [^ T(x, y, z, v) dtd. 


p(z) dz p(z)oz 


Int J Adv Appl Sci, Vol. 11, No. 4, December 2022: 345-355 


Int J Adv Appl Sci ISSN: 2252-8814 o 353 


Approximations of the second and higher orders of components of displacement vector could be 
determined by using standard procedure. The equations for the required components after the standard 
substitution takes the following form: 


a? C 3 5E(z) a? (x,y,Z,t) E(z) 07 uz y(x,y,z,t) 
p(z) oe = {K(z )+ z jJ ux Goy,z + ike - Zz. j 1y + 


6[1+0(z) 0x? 3[1+0(z) axdy 
87 ury(x,y,2,t) E(z) [ene 07u47(x~,y,z,t) . OTGoyzt) x 
ðxðy 2[1+0(z)] 0y? 0z? Ox 
E(z) } Pian 
P KOP) + KO + 3[1+90(z)] OxOz 
( ) 0?u;y (xy z,t) _ E(2Z) [Ease 07u4x(x,y,z,t) _ OTGoyzt) 
p 8t? ^ 2[1+0(z)] ax2 axdy ay 
E(z) [eee Ou4zGQoy,z.t) ] 0^uiy yz t) 
% K(z)B(@) tp Oz iS Oz T oy + 0y? ^ 
5E(z) } ( || E(2) uen aury GG z,t) 
% A T K(z) T K(@) 6[1+0(z)] ðyðz + K(z) ðxðy 
(z ) aso z o ..  E(z) [eene O^ujz(xy,zt) , O^?u4x Gy zt) + 
PY 0t? ^ 2[14o(2)] 0x? 0y? Oxdz 
ugs d 4 2 [kc ) [pasen ds Muse 4 Tee + 
E(z) 9 | Eom t) E Ou4xCoy,zt) = Oury(x,y,Z,t) z Ou4z oy. z,t) = 
6[1+0(z)] dz Oz Ox oy Oz 
2 Ou4x Gy,z,t) = OuiyG y zt) " Ou4zCGoy,z.t) ] E(z) i K(z JLE je usn rene 
Ox oy Oz 1+0(z) 


Integration of the left and right sides of the above relations on time f leads to the following result: 


5E(z) 1 
uyr) = = (KG + ONS h fp tas Go z dede +- {K (2) 
E(z) E(z) 
"spes EN. h Ury x, y, Z, T)dtdd + 2009 se f: Uy x, y, z, T)dtdð + 


+f K uz(x,y,Z, z)drdó| — f ig u4 (x, y, Z, T)dtdd (K (z)+ 


1+0(z) Hs Oxdz 


E(z) B8 
xe K(z 2) ala (s T (x, y, Z,t)dtdd -2p f? ux x, y, Zz, T)dtdd x 
5E(z) E(z) 
ASK (z ) E dies - {K (Z ) z 3[1+0(z)] JA R Ury (x, »2 t)dtdd x 
A E(z) 


r ETT ET T b ds Uy (%, y, z, T)dtdd + — E 1s uz yz z)drd$| - 
make ) +2 o Ep Jy tz% y, z, dtd + uox + K(z) 2 x 


3[1+0(2)] 

ee f? T (x, y, z, tT)dtdd 
u5y (x, y,z, t) = E ETIN Urix (x,y,Z,T)dtdd + ah D ux yz, z)d«dé| + 

Eg E 2i i ta Ne) dede p) eso t kc) A 

Sul IM ux Y, Z, T)dtdd + S x ates lee s u4y x, y, z, T)dtdd + 

e P e yz t)drdð |) -KEE h h Tœ yz dedo - {Oo - 
-KO aaah j uy G5, y, z, t)dtd8 — e [zc js ux x, y, Zz, T)dtdd + 

zah l uix (X,Y,Z, z)drdó| — E. z im K(z rons ih T (x,y,z, T)dtdd — Eo 

Xy a u, (x, y, Z, T)dtdd — EDT Ta Ds Uy x, y, Z, T)dtdd eus t 
TU) ren E hs fs Ury(x, y,z, Tdtdd tL is Uz (X,Y, Z, x)d«d$||] x 


1 E(z) 
x 2p) -fk (z ) B 6[1+0(z)] j pah k Ury (x, WZ, idi + Moy 


à 
uz(x, y,zt)- x ex f U4, (x, y, z, t)drdO +5 a T uiz (x, y, Z,T)dtdd + 
ð? o p9 1 
Feo I ux x, y, Z, t)dtad + TM IM uiy Gy; Z, mE uz) + a x 
ð ð o pV ð ro rpt 
x te ) Es h e AA Que t 21, Sy x05 y, z t)dtdd + 
(z) ð po po 
IN Tus OZ, t)drdð|} + E — 2l a [e 7 f; Sy Uz y, z 0)d1d8 — 


1+0(z) 
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ð po p ð po p ð po pV 
E So Uy x(x, y, Z, T)dtdd c J Ury x, y, z, T)dtdd Er So us2(x,y, 2, t)drdo]} - 


Bz) ð po p8 
—K(z) seines J, TG yz, 1) drdO + uoz- 


In this paper, we calculate concentrations of a mixture of gases and a diffusing material in the 
considered heterostructure, distribution of temperature, and components of displacement vector as the 
second-order approximations in the framework of the method of averaging functional corrections. This 
approximation is usually enough good for obtaining qualitative conclusions and obtaining some quantitative 
results. The results of analytical calculations were verified by comparing them with the results of numerical 
simulation. 


3. DISCUSSION 

In this section, we analyzed the diffusion mixing of the heterostructure layers during annealing with 
account for the relaxation of mismatch-induced stresses. Typical distributions of diffusant concentrations in 
the considered heterostructure are shown in Figure 2 for the different continuance of annealing time. The 
conclusions, in this case, are standard: an increase in the duration of annealing leads to a more homogenous 
distribution of the diffusant. At the same time mixing of heterostructure materials leads to decreasing in 
mechanical stresses as shown in Figure 3. It should be noted that the porosity of the buffer layer leads to 
increasing this effect. 


1.0 4 
0.8 4 
1 4 
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© 2 
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3 o2 - 
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0 Z L. 0.0 a 
r4 
Figure 2. Spatial distributions of diffusant in the Figure 3. Normalized dependences of the component 
considered heterostructure at different values of the of the displacement vector uz on the z coordinate (a is 
continuance of annealing time. Increasing of the thickness of the buffer layer). Increasing of 
number of curves corresponds to increasing the number of curves corresponds to increasing the 
continuance of annealing time continuance of annealing time 


4. CONCLUSION 

In this paper, we analyzed the effect of processing the substrate, which precedes the growth of new 
epitaxial layers, on the properties of the grown heterostructure. It has been shown, that growth of the new 
epitaxial layer on the buffer layer after preliminary (before the start of growth) annealing decreases the value 
of mismatch-induced stresses. An analytical approach has been introduced for the analysis of mass and heat 
transfer in a multilayer structure with account mismatch-induced stress. 
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